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We analyze interacting one-dimensional bosons in the continuum, subject to a periodic sinusoidal potential 
of arbitrary depth. Variation of the lattice depth tunes the system from the Bose-Hubbard limit for deep lattices, 
through the sine-Gordon regime of weak lattices, to the complete absence of a lattice. Using the Bose-Fermi 
mapping between strongly interacting bosons and weakly interacting fermions, we derive the phase diagram in 
the parameter space of lattice depth and chemical potential. This extends previous knowledge from tight-binding 
(Bose-Hubbard) studies in a new direction which is important because the lattice depth is a readily adjustable 
experimental parameter. Several other results (equations of state, energy gaps, profiles in harmonic trap) are pre- 
sented as corollaries to the physics contained in this phase diagram. Generically, both incompressible (gapped) 
and compressible phases coexist in a trap; this has implications for experimental measurements. 



Introduction — Ultracold atomic gases provide us with 
unprecedented experimental opportunities for creating inter- 
acting quantum many-particle systems [1, 2|, both in lattice 
and in continuum situations. In response, theoretical lattice 
and continuum models have been widely analyzed in the past 
decade and half. However, it is experimentally straightfor- 
ward to interpolate between the two cases, through periodic 
optical potentials whose amplitude can be made to vary from 
zero to a large value |3 |. The heating or cooling due to ramp- 
ing up the strength of the periodic potential is well-studied 
However, the potential well depth is generally not considered 
as a parameter tuning the phase diagram, except indirectly via 
effective Hubbard-model parameters which are strictly valid 
only in the deep-potential limit. The study of complex many- 
particle systems as a function of the periodic potential depth, 
interpolating between continuum and tight-binding limits, is 
a direction with many interesting effects whose exploration 
is just beginning. The measurements of Ref. |3| are promis- 
ing first steps in this direction. This experiment focuses on 
one-dimensional (ID) physics, where dimensional confine- 
ment enhances correlation and interaction effects. 

In this work, we address the system implemented in the 
experiment of Ref. 13], namely, one-dimensional strongly in- 
teracting bosons in the presence of a variable-strength opti- 
cal potential. We provide the most fundamental information 
needed for understanding the physics in this system, namely, 
the phase diagram. 

Like many other systems, ID bosons have been widely 
studied in the tight-binding and continuum limits |5|, but not 
much for intermediate depths. The deep-well limit is de- 
scribed by the Bose-Hubbard model, which is a tight-binding 
model in the sense of each well being described by a single 
mode. The ID Bose-Hubbard phase diagram is well-studied 
theoretically |l5Hll> ™d the model has been realized often 
experimentally |[I]|2]|5]. In the other limit, i.e., without a 
periodic potential, the system is the continuum Lieb-Liniger 
gas [9 1, described by the Hamiltonian H ~ /dxj'^ + 
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FIG. 1: (a-c) Schematics of the /i-7 phase diagram for weak, in- 
termediate and deep lattices. The Mott regions are shaded. The 
deepest lattice is taken to be described by a Bose-Hubbard model 
|7|. (d,e) Calculated \i-s phase diagram for the strongly interacting 
(7~^ = 0) limit. Each Mott region corresponds to a different inte- 
ger filling V. The dashed blue line is the energy of a single particle 
placed in an empty system, i.e., the chemical potential for v = 
(f,g) Inferred schematics for finite 7~^: phase diagram and critical 
depths, (h) Structure of the 3D phase diagram (in /i-7-s space) show- 
ing the 'Bose-Hubbard' region (previously studied) and the large-7 
plane (determined in this work). 



gJ2i<j ~ ^j)- Th^ dimensionless interaction parame- 
ter is 7 = gm/ [h^n), with m the mass of bosons and n the 
one-dimensional density. We are interested in large 7, in or 
near the so-called Tonks-Girardeau (TG) regime. The TG gas 
has been intensively studied theoretically ll5l iT0H23]| . Some 
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calculations exist in the weak lattice limit, via mapping to a 
sine-Gordon field theory Il24lj26l . The TG gas has also been 
experimentally realized in a lattice version with low filling 
127 J . Ref. [3J studies strongly interacting ID bosons in the 
continuum, with the addition of a periodic potential of vari- 
able strength. Motivated by this, as well as by the dearth of 
theoretical studies of the effect of the potential depth, we ana- 
lyze a Tonks-Girardeau gas in a periodic potential of variable 
depth. The periodic potential is V{x) = sEi^ sm^{qx), with 
lattice period n/q. The parameter s gives the potential depth 
in units of the recoil energy, Eji — . 

We use the Bose-Fermi transformation (T0 \ to map strongly 
interacting ID bosons (7 — > 00) to free fermions. The prob- 
lem with a periodic potential then turns into a band-structure 
problem, which we address using the Bloch theorem. This 
rather simple setup allows us to extract a remarkable amount 
of information. The most prominent result is the phase dia- 
gram in the chemical potential (/i) versus well depth (s) plane, 
i.e., we map out parameter regions which are Mott phases and 
those which are superfluid phases [Fig.[TJd,e)]. This phase di- 
agram, whose structure to the best of our knowledge was not 
previously known, is exact and detailed in the 7 — > 00 limit, 
and also provides a significant amount of information about 
the phase diagram for finite but large 7. 

We also present excitation gaps in the different Mott phases. 
Defining a "filling" variable as the average density in a well, 
we present equations of state, i.e, /i(i^) curves, for various well 
depths. In addition, we calculate the effect of an arbitrary- 
strength periodic potential in the presence of an overall har- 
monic trap. This is important because, like Ref |3|, fore- 
seeable experiments are likely to be performed in harmonic 
confinement. Going beyond the infinitely interacting case, we 
also provide results for large but finite interactions. 

The phase diagram — Fig. [TJa-c) illustrate the phase di- 
agram in the plane of chemical potential /i and (inverse) in- 
teraction strength, for weak, intermediate and strong lattice 
strength s. The deep-potential limit has been represented as 
the well-studied Bose-Hubbard model (top right). Away from 
this limit, the phase diagrams are expected to be of the forms 
sketched, but are not presently known in detail. It should be 
possible to determine parts of the small-s phase diagrams us- 
ing the sine-Gordon model, but full results do not yet exist in 
the literature. 

Our calculation determines the small 7"^ limit, i.e., the 
region touching the vertical axes in each of the top-panel 
schematics, for all values of the lattice depth s. This can be 
put into context in terms ofthe 3D phase diagram of Fig.[TJh): 
our calculation maps out a new plane in pi-'y-s space. 

For 7^^ =0, the bosonic phase diagram is obtained from 
the free-fermion band structure. A filled fermionic band cor- 
responds to integer filling in the bosonic system. Thus, the 
fermionic band-gap regions correspond to Mott phase regions 
in our system. The resulting exact phase diagram is shown in 
Fig.[TJd,e). This phase diagram is our main result. 

The s — > cx) region satisfactorily reproduces known fea- 




FIG. 2: Left: Equation of state for various optical lattice depths. 
Solid curves correspond to 7 — )■ 00, dashed curves to 7 = 11. Right 
top: energy gaps Eg m 1^ = 1, 1/ = 2 and 1/ = 3 Mott phases (top 
black, middle blue and bottom red) as a function of lattice depth. 
The gaps scale as s" at small s and as s^^^ at large s. Right bottom: 
calculated 1/ — 1 energy gap (line) shown with experimental values 
from trapped-system measurements (3|- 



tures of the Bose-Hubbard limit — the Mott regions touch, 
killing superfluid regions, and each occupies an equal range 
of /i. In the shallow "sine-Gordon" limit (s — > 0), the Mott 
regions shrink to tiny slivers. For 7^^ = these Mott sliv- 
ers continue down to s = 0, so that an infinitesimal periodic 
potential can pin a gas if the filling allows commensurate pin- 
ning. In the s — limit, the successive Mott slivers appear at 
linearly increasing gaps of ^ from each other. 

In Fig. [TJf) we sketch the expected phase diagram for large 
but finite 7. The Mott regions now start at some finite crit- 
ical lattice depth Sciit- The Mott regions are expected to be 
less robust for larger fillings because of the smaller gaps and 
stronger perturbative corrections at larger 1/ (see below). It 
remains an open problem to determine the Scrit (7) curves for 
1/ > 1, sketched in[TJg)- The i' — 1 curve can be obtained 
from a sine-Gordon description lO. The other curves should 
be be progressively higher [Fig. [ijg)] but are otherwise un- 
known. 

Bose-Fermi mapping — To obtain the wavefunction of 
either N hardcore bosons or N noninteracting fermions 
in an arbitrary single-particle ID potential, we take 
the first N single-particle eigenfunctions Ui{x). The 
many-body wavefunction is then '(/;/(a;i, 0:2, x^v) = 
(— 1)^ Wi(Pxi) (where P denotes a permutation) 
for fermions and 7/'f,(a;i, X2, xjv) = X]p Si "i(-P^i) 
bosons |10|. Thus, we need to solve the single-particle 
Schrodinger equation with potential V{x) = (a;^/2£jj.) + 
sEjism'^{qx), where £tr is the trap length for the harmonic 
trapping potential. For a uniform (non-trapped) system, the 
solution is known to be expressible in in terms of Mathieu 
functions |28,29|. 

Once the lowest N single-particle orbitals Ui{x) have been 
calculated, the one-dimensional density can be shown to be 

n{y) = / dx2-..dxN\tpiy,X2,...,XN)f ^^uf{y) . (1) 
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The "filling factor", the number of particles divided by the 
number of wells, is (for uniform systems) the integrated n{x) 

over one period: z/ — Jy'^^^'^ dxn{x). For non-uniform sys- 
tems we assign to each well j the filling ^{j), defined as the 
integral of n{x) within that well. 

Energy scales — For deep wells (large s), expanding 
sain^lqx) w sq^x^ shows that the lower energy levels have 
equal spacing cx y/s. This implies equispaced bands, which 
explains the equal widths (cx y/s) of the Mott regions. On 
the other hand, for s — > 0, zero-point energy arises due to 
confinement by the neighboring particles. At filling ly this dis- 
tance is Tr/i^q and therefore the relevant energy is v^Ej^. This 
explains the positions of the Mott slivers for s — and the 
linearly increasing distance between successive Mott regions 
[Fig.[3d,e)]. 

Perturbation theory away from the Tonks-Girardeau point 
— Near the TG limit, i.e., for finite but large interactions, the 
ID Lieb-Liniger gas can still be mapped onto a ID weakly 
interacting fermionic problem. The fermionic interaction is 
unfortunately not simple. We use the form lfT2HT6l 
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S"{x) = -gfS"{x). 



(2) 



This is to be interpreted as a perturbation to the fermionic 
state. Thus the energy shift for the bosonic ground state is 
SE = {ipf \ V \ipf)- After some manipulation we obtain 
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with (72'" denoting the m-th (n-th) derivatives with respect 
to the first (second) argument of the pair correlation function 
g2{x,y). For the TG gas, g2{x,y) = n{x)n{y) - |A(x,?;)|^ 
where A(a;,y) = 'Yl,i'^i{x)ui{y), the sum running over the 
first N eigenstates EOl . 

Equation of state — Fig. 2 (left) shows the chemical poten- 
tial versus filling, for several different depths. We call these 
jjL versus v curves the "equations of state". As the ground 
state is constructed by successively occupying single-particle 
states, the ij.{v) curves are given by the single-particle energy 
dispersion curves (energy e versus momentum k), with the in- 
terpretation V = k/q. The Mott regions are incompressible 
in the sense that the filling does not change with chemical po- 
tential, hence there are segments of these curves which are 
vertical jumps at integer values of v. The topography of the 
phase diagram — thin Mott regions for small s, thick Mott 
regions for large s — is visible in the structure of these ^(i^) 
curves. 

We also show the perturbative correction to the equation 
of state curves for 7 large but finite. Specifically, we display 
7 = 11, corresponding to some of the experimental data in 
Ref. 131- For small and moderate lattice strengths (,s < 10), 
the corrections are minute. For s ^ 10, perturbative results 
are difficult to interpret and thus not displayed. 



Excitation gap in Mott phase — Fig. 2 (right top) shows 
the first few energy gaps as a function of the trap depth when 
the system is in the Mott phase. The gaps correspond to the 
vertical jumps in the /i(i') curves (Fig.l2]left) or, equivalently, 
to the width of the Mott regions in Figll] We find Eg/En — 
) ^ '^I'l^f)!' with and bi, the characteristic values 
of Mathieu functions |29|. For small s, the v-th gap grows as 
s'^ , e.g., as s/2 and as s^/32. For v ~ 1 this is in agreement 
with sine-Gordon results Il3l l24l|251 . For s 3> 1, all gaps cross 
over to ~ s^/^, in agreement with the energy-scale arguments 
given previously and consistently with Fig.[TJe). 

The v ^ 1 gap has been experimentally measured in lO 
through modulation spectroscopy. In Fig.[2](right,bottom) we 
show the experimental data values (for 7 = 11) together with 
our exact calculation for 7 = 00. The perturbative result for 
7 = 11 is indistinguishable from the TG line. Our results for 
trapped systems, below, show that a uniform-system gap can- 
not be expected to coincide with measurements on a trapped 
system with significant inhomogeneity. 

An overall harmonic trap — In Fig.|3]we show the exact 
density profiles of a harmonically trapped Tonks-Girardeau 
gas in an optical lattice, for various lattice depths. We show 
both the densities n{x) as a function of space and the fillings 
v{j) as a function of well index. The fillings v are obtained 
by averaging n{x) over each well. 

The fillings (bottom row) show a wedding-cake structure 
familiar from the literature on Bose-Hubbard physics in traps. 
It is noteworthy that we have obtained this from averaging 
densities in a continuum model, and not from a tight-binding 
model like the Bose-Hubbard. The v{i) curves can be under- 
stood from "local density approximation" arguments based on 
the phase diagram presented in Fig.[T] Going from deeper to 
shallower modulations, the Mott regions grow thinner [Fig. 
[TJd,e)], and correspondingly there are smaller plateau regions 
in the density profiles. This resembles the effect of going from 
stronger to weaker interactions in the Bose-Hubbard model. 

The continuum density profiles, n{x), are less familiar from 
other contexts. The density oscillates at the lengthscale of the 
well size, so we have highlighted some features in insets. For 
'fillings' larger than one, n{x) within each well typically has 
multiple peaks. The height of the peaks in the v = 2 plateau 
region is larger than, but not twice, that in the v = \ region. 
Only the filling v{j), i.e., the integrated density, has integer 
values at the plateaus. Another feature is that the densities do 
not completely vanish between the wells, especially at small 
s, in neither superfluid nor Mott region. 

In Fig.|4] we show profiles tailored to the experimental sit- 
uation of Ref 13 !, where the central tubes are reported to have 
about 60 bosons with maximal filling around v « 1.2. In 
a non-trapped (uniform) system at integer v, an arbitrarily 
weak modulation potential pins the whole system in a Mott 
state. Our results in Fig. |4] show that the situation is more 
complicated in a harmonic trap; v is now strongly position- 
dependent. In particular, in the s < 1 regime, the part of the 
cloud which is 'pinned' at = 1 is a small fraction of the to- 
tal. Therefore, identification of modulation spectroscopy mea- 
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FIG. 3: (Color online.) A harmonic trap in addition to the modulated potential. Top panels: densities n{x), as a function of continuous 
position X. Bottom panels: fillings as a function of discrete lattice index j. Blue curves are for 181 bosons. Lower black curves, where 
present, correspond to 31 bosons. The trap length is SO/vr times the modulation period. Insets highlight features discussed in the text. 
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FIG. 4: Left: Density profiles in a trap, with parameters correspond- 
ing to the central tube of the experiment of Ref. |3|. Right: Molt 
fraction, defined as the number of sites with — 1| < e divided 
by the total number of bosons, with e = 0.02(0.05) for the lower 
(upper) curve. 



surements lO in this regime with a "Mott gap" is nontrivial. 
Further investigation incorporating inhomogeneity consider- 
ations is therefore called for. While Fig.j4] shows 7 = oo 
results, the perturbative calculations in Fig.l2pndicate that cor- 
rections to the profiles will be tiny for interactions as low as 
7 ~ 10. 

Discussion; Open questions — We have studied strongly 
interacting bosons at and near the hardcore limit for periodic 
potentials of arbitrary depth. We have presented the phase di- 
agram of the system in the jjL-s plane, and quantitatively char- 
acterized various related aspects. By providing an exact and 
detailed map of one plane of the ji-^-s space [Fig.[TJh)], we 
give a good first idea of the topography of this 3D phase dia- 
gram, and pave the way for future explorations of the complete 
3D phase space. 

Not surprisingly, the Mott regions are less dominant for 
weaker lattice strengths. This has clear consequences for 
trapped systems — the fraction of the system locked at in- 
teger values is small at small s. This may explain why the 
experimentally measured values for the v — \ gap fall below 
the theory for the uniform system. Heuristically, a significant 
part of the system is compressible (gapless), which may be 



expected to reduce the measured gap. Since the modulation 
spectroscopy employed in Ref. fS\ involves complicated tem- 
poral dynamics, this highlights the need for studies of time 
dependent properties of the non-uniform system with a weak 
lattice. 

Other open questions raised by this work include the ef- 
fects of non-infinite 7, beyond the perturbative calculations 
presented here. For s ^ 10, we find that the perturbative 
corrections at 7 ss 10 lead to the fermionic dispersion rela- 
tion becoming non-monotonic, so that the Bose-Fermi map- 
ping would involve careful (re-)interpretation. At small s, the 
perturbative calculation could in principle give estimates for 
the Scrit value where the Mott slivers end. (For 7^^ = the 
Mott slivers extend to s = 0.) In practice, we have found that 
the perturbative corrections to the gap are so small that exten- 
sive high-precision computations would be required for such 
a calculation. A Monte Carlo or DMRG study, analogous to 
the 3D study of |30|, would be useful in this regard. 
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